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PREFACE 


aA 
bee 


This paper presents the results of a study of Constant 
Absolute Vorticity Trajectories applied to isentropic surfaces. 
The objectives of this study were: first, to show that Constant 
Absolute Vorticity Trajectories are theoretically better applied 
to forecasting future positions of parcels on isentropic surfaces 
than on constant pressure surfaces as they are now applied; sec- 
ondly, to show quantitatively the actual improvement of the tra- 
jectories on isentropic surfaces by comparing statistically 
results of actual forecasts from both types of charts; thirdly, 
to suggest how the improved technique can aid forecasting, par- 
ticularily upper air prognosis. 

This work was undertaken as the thesis requirement for the 
degree of Master of Science in Aerology at the U. S. iaval Post- 
graduate School, Monterey, California, during the acadeniic 
year 1952-1953. 

The author is indebted to Associate Professor George J. 
“altiner of the Aerology Department for his very valuable assis- 
tance and constructive criticism during the investigation. Ue 
also wishes to acknowledge the assistance of Professor A. Boyd 
Mewborn of the Mathematics Department and the author's wife, 
Pauline, who contributed so much help with the laborious task 


of decoding, computing, and entering isentropic charts. 
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I. IPTRODUCTION 


The meteorologist is often called upon to forecast for periods 
in excess of twenty-four hours for which simple extrapolation is 
generally insufficient. Gne method of forecasting entails prepa- 
ration of prognostic surface charts and preparation of forecasts 
from these charts. As most meteorlogists Imow, making accurate 
prognostic surface charts is easier said than done. One of the 
important tools used in estimating the position and intensity of 
surface systems is a prognostic 500 millibar chart. At present, 
much attention is being paid to the problem of upper air prognosis, 
and considerable information has been published concerning fore- 
casting the 500 millibar surface, the most recent being Forecasting 
in the Middle Latitudes by Richl and collaborators. It is hoped 
that this paper may shed some light on a facet of this problem and 
may contribute to better prognosis and forecasting. 

One aid to prognosis of upper air charts is the construction 
of forecast air parcel trajectories based on the principle of con- 
servatism of the vertical component of absolute vorticity ,(¥ +f ) F 
where ay is the vertical component of relative vorticity and + 
is the vertical component of the coriolis force known as the 
coriolis parameter. tossby {e} initiated the vorticity concept 
by showing that under certain restrictive assumptions, the follow- 


ing relationship holds: 


(f+f) = constant. 


(1) 
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The assumptions referred to above are as follows: 

a. The atmosphere is barotropic. 

b. The atmosphere is a homogeneous, incompressible fluid. 

c. Motion is purely horizontal. 

d. Friction forces are neglected. 

e. There is no horizontal divergence. 

Rossby fo, pp 268-2897 then integrated his fundamental equation 
and was able to obtain trajectories for particles which conserve 
their vertical component of absolute vorticity. This trajectory 
is known as the Constant Vertical Component of Absolute Vorticity 
Trajectory and will be referred to as a CAV trajectory. 

The assumptions necessary to derive the vorticity equation 
according to Rossby are obviously quite restrictive. Starr, as 
editor of the Journal of Meteorology (June, 1945), suggested that 
CAV trajectories mignt better be depicted on an isentropic chart 
than on a constant level chart as there are no solenoids on an 
isentropic chart. Namias [s, pp 372-374] states that from obser- 
vations available isentropic surfaces appear, to the first approxi- 
mation, to be substantial surfaces, i.e., surfaces that contain the 
same air particles from day to day. Influences of non-adiabatic 
nature are appreciable over long periods of time, but are usually 
insufficient to disrupt the fundamental isentropic flow from one 
day to the next. The derivation of the vorticity equation for flow 
on an isentropic surface requires fewer assumptions. Such restrictions 


as purely horizontal flow, a barotropic atmosphere, and a homogeneous 


(2) 
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incompressible atmosphere may be removed. Now, an isentropic 
surface is a surface of constant potential temperature. In 
meteorlogical processes, potential temperature is one of the more 
conservative elements by which parcels of air may be identified 
from time to time. Since the atmosphere is normally stable, 
potential temperature increases with height and the atmosphere 
ray be considered to consist of an infinite number of isentropic 
surfaces. Sir Napier Shaw fro, p 263} first suggested that charts 
of isentropic surfaces should be drawn as the motion of air is best 
resolved on isentropic surfaces. It will be shown in this paper 
that the prognostic trajectory based on the conservatism of the 
vertical component of absolute vorticity is theoretically and 
practically better suited to use on isentropic charts than on 


constant pressure charts as it is now applied. 


(3) 
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Ti. THEORETICAL INVESTIGATION 


The motion of particles on a substantial surface can be 
described with Lagrangian methods. The form of the vorticity 
equation for flow on an isentropic surface derived here using the 
Lagrangian method of attack is due to Yaltiner”. Let x,y, Z andt 
be normal Cartesian coordinates and lec EY, 7 @a Y be 


Lagransian coordinates where % and Y are coordinates of the 


horizontal projection of the isentropic surface, 7 = 2 andy 


refers to time for the isentropic surface. Then 
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It follows from the above equations that 


} 
See a (2) 
fos oY 


the stream function for an isentropic surface. 


where Fuel t+qe ; 


In a sinilar way we find 


gu . Oe dude , Ov — dv _Qv Oz 
Oy OOK 25K Ox OX 92 dX (3) 
Ou = Du ~ du Iz 9 oe Ow _ dv dt (4) 


*Associate Professor C. J. Haltiner, U. S. Naval Postgraduate School 


(4) 
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Further, since the motion is assumed to 


ue = S& 4udE +v OZ. 
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Substituting in the equations of motion 
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the expressions 1, 2, 3, 4, and 5, we obtain 


du pudu vow — ~_OF + fv 
Yi) 
as Hage ~ Fw {fo\ 


Differentiating (7) with respect to Y and (8) with respect to X 


and subtracting, we obtain 
[5 Ow Jv tie Qe pd + Ow de pu d wu _oe ov 
oy¥y/OLXY KX OY OX OFYOY oeY of aY*t OX oY 


oy ak dv 2 of OF 
—-HASM -vew -— Sore => Sviet — —— 
om” YOK OM ay oY a Ox a Sy OX 


Collecting terms 
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and introducing vector notation we obtain 


~ (SLL) +V%K)-L Be VaV EF + ABV 42 
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where the operator Vo= Ho 71 2 


Finally, we may write this as 


LK +f +(2+4)Yov =o (9) 


Since p~p=t for any particle. Now, assume the divergence term Wj,V=0 


Then we may write the vorticity equation in the familiar form 


AL? a =O or ( f + ) = comsfanl. (10) 


While this form is similar to the form derived by Rossby, there 
are some fundamental differences. Parcel motion is not limited 
to the horizontal, but is three dimensional. Moreover, a barotropic 


and incompressible atmosphere is not required. 


(6) 





IIL. TECHNIQUES OF INVESTIGATION 


The relative vorticity may be written in a form similar to 
the usual expression for the vertical component of vorticity in 


natural coordinates 


— Vv. by 
k= ON 


where R 5 and N are measured on the isentropic map. If we select 
a parcel in a broad uniform flow or in the axis of a jet stream, br =o) 
Further, if the parcel is at an inflection point in the flow, M-o. 
Under these conditions the CAV trajectory can be computed. It ; 
will be noted that R s ‘defined above is not exactly equal to the 
projected radius of curvature of the streamline on the isentropic 
surface. However, the approximation involved here is normally 

less than those made in the actual computation of CAV trajectories 
as described by Fultz [4 p i} and revised by Wobus [12] - While 
approximations have to be made, it must be remembered that the true 
worth of a technique lies in the statistically proven results 
obtained. 

The simplest method of obtaining a CAV trajectory is through 
the use of a mechanical device such as the so-called "Wiggle-Wagon" 
devised by Wobus [ 12 ] and currently in use in the WBAN Analysis 
Center, Washington, D. C. Original studies of the CAV trajectory 
assumed a plane earth. This assumption entailed errors due to 
distortion when compared to trajectories along a spherical earth. 


The Wobus "Wiggle-Wagon" corrects this defect. To aid meteorolo- 


gists who do not have access to the "Wiggle-Wagon", the U. S. Navy 


(7) 





Bureau of Aeronautics Project ARowA [11] has published a table 
for computing CAV trajectories made from trajectories traced by 
the "Wiggle-Wagon". This table will be used to compute CAV tra- 
jectories in this paper. For application, it is necessary to know 
speed, direction, and the latitude of the parcel, all at an inflection 
point in the flow. 

The next step in the investigation was to compute trajectories 
on isentropic and 500 millibar charts and test the results statisti- 
cally. A rather complete compilation of upper air and surface data 
are available in the Northern Nemisphere Surface and 500 liillibar 
Charts Series published by the U. S. Weather Bureau. These remarkable 
publications not only have a daily series of logically analysed surface 
and 500 millibar charts, but also contain a compilation of all data 
used to plot these charts. Data were taken from the period of 
January-February 1949. From this data, a set of thirty isentropic 
charts was plotted. The 303K isentropic surface was selected because: 
(1) this surface was often close to the so-called "level of non- 
divergence’! described by Bjerines [2] ; (2) this surface often con- 
tained much of the 500 millibar wind speed maximum axis (jet); 
(3) the 303K isentropic surface was usually sufficiently far fron 
the surface to remove effects of surface friction and turbulence. 
It was desired to select an isentropic surface containing much of 
the 500 millibar isotach maximuni axis so computation of CAV tra- 


jectories could aid in the prognosis of the 500 millibar surface. 


(8) 
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Inasruch as the derivation of CAV trajectories on an isentropic 
surface assumes the divergence term to be zero, the selection of 

an isentropic surface close to the level of no horizontal divergence 
may be desirable. The effect of the divergence term is discussed in 
Chapter IV. 

After the heights of the 303i! isentropic surface were plotted, 
all available wind data pertinent to the surface were plotted. Where 
Wind observations were available from stations other than radiosonde 
stations, an interpolated value of the heignt of the isentropic surface 
and the wind nearest that height were entered. Because of the fairly 
dense network of upper air wird observing stations in the United 
States, sufficiently accurate streamline analysis was possible with- 
out the time-consuming computations of stream function values. 

Streamline analyses were made for 0300Z each day in the period 
from 6 January, 1949, to 22 January, 1949, and from 10 February, 1949, 
to 20 February, 1949. For each chart, a number of CAV trajectories 
were conputed. The following conditions had to be met before choosing 
the initial point: (1) the point had to be in or very near an inflection 
point in the flow (undergoing no effects of curvature); (2) the point 
had to be in or near the axis of the isotach maximum (jet stream); 

(3) the point must be at an actual point of wind observation to obtain 
correct initial wind speed and direction. This procedure netted from 
one to four points suitable for computation of CAV trajectories per 


chart. 


(9 
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For comparison with each CAV trajectory plotted on the 
isentropic chart, a companion CAV trajectory was computed on the 
500 millibar chart of the same time. It was not too difficult 
to find a companion inflection point on the 500 millibar chart as 
the isentropic surface was usually close to the 500 mnillibar surface. 
After the CAY trajectories were constructed on each chart, forecast 
wind velocities were determined by moving the parcel along the tra- 
jectory a distance equal to the initial wind speed multiplied by 
the number of hours to forecast time. The forecast times selected 
were twenty-four, forty-eight, and seventy-two hours. These fore- 
cast winds were then plotted on the verifying charts and verified. 
On the isentropic chart, verification depended on the forecast 
position falling within the observational network. If the forecast 
position fell within the network, but not on a particular station, 
the forecast winds were verified by linear interpolation between 
nearby stations. However, care was taken to use only stations 
where linear interpolation would be reasonable. On the 500 millibar 
chart, forecast winds were verified in a similar manner if the fore- 
cast position fell within the observational network or on a reporting 
Station. If the forecast position on the 500 millibar chart fell on 
neither of these places, TeeconSeune were assumed to be streamlines, 
and gradient winds were measured to verify. Since verification on 
isentropic charts depended on forecast positions falling in the North 


American network or on a particular station, the number of possible 


(10) 





verifications fell off drastically after forty-eight hours. lHowever, 
with hemispherical 500 millibar charts available, forecasts made on 
the 500 millibar chart could always be verified. All 500 millibar 
CAV trajectories were carried out to seventy-two hours. 

In verifying forecast winds, observed wind direction error 
was arbitrarily selected negative if the observed wind direction was 
to the right (clockwise) of the forecast wind direction and positive 
if the observed wind direction was to the left (counter-—clockwise) 
of the forecast wind direction. Similarily, wind speed error was 
Selected negative if observed wind speed was less than forecast and 
positive if the wind speed was greater than forecast. All forecast 
winds and verifications were tabulated in a form suitable for 
statistical testing. 

Histograms of wind speed error and wind direction error were 
prepared for twenty-four and forty-eight hour forecast times for both 
the isentropic and 500 millibar charts. These histograms are depicted 
as Figures 1 through 4. Wind speed errors were grouped into cells 
each having a cell interval of 10 knots error. The cell midpoints 
are fort O, £190,420 etc., knots error. Since the cell boundaries 
are exactly halfway between the cell midpoints, some of the obser- 
vations fell on the cell boundaries. in this case, it was deemed 
best to divide such observations between each adjacent cell. Less 
trouble was encountered in grouping wind direction error. All wind 
direction errors were in increments of 10 degrees and values of 0,+/0 


+20, etc. degrees of error were taken as cell midpoints. 


(a) 
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After the data were grouped in histograms, normal curves were 
fitted in accordance with the method described by Hoel [s, pp 191-194]. 
After obtaining these fitted normal curves, XK tests were performed to 
determine if each sample of data could reasonably be from a normal 
population. In all cases the sample distributions could reasonably 
be said to be from normal populations, although in two cases the 
value of e was close to the critical value. Restrictions on the 
test were met by grouping cells with too few frequencies together. 
After showing normality of the data, the distributions of data from 
the isentropic charts were compared to similar distributions from 
the 500 millibar charts by the use of the "F" test. The "IF" test 
was made to determine if compared data distributions from the 
isentropic and 500 millibar charts could statistically be from the 
Same parent normal population. The purpose here is to show that a 
given distribution of data from the isentropic charts is not from the 
same normal population as the similar data distribution from the 500 
millibar charts. If there is a significant difference, then inspection 
and intuition will reveal which of the two types of charts are better 


suited to forecasting by means of CAV trajectories. 


(12) 





IV. “RESULTS 


It is well to look first at the data in "everyday" terms to 
find what kind of results a forecaster may expect. If it is arbi- 
trarily decided that a wind direction forecast will verify as a "hit" 
if it is within 20 degrees of observed direction, then we find that 
96.4% of the twenty-four hour wind direction forecasts from CAV tra- 
jectories on isentropic charts verified as “hits'. By comparison, 
76.35 of similar forecasts from 500 millibar charts verified as "hits'', 
With forty-eight hour wind direction forecasts, 81% from isentropic 
charts verified as "hits" and 41% from 500 nillibar charts. 

If it is arbitrarily decided that 15 knots of wind speed error 
is not too much to verify as a "hit", we find that 763 of twenty-four 
hour wind speed forecasts fron trajectories on isentropic charts 
verified as “hits”, Of the data from the 500 millibar charts, 42.4. 
of twenty-four hour wind speed forecasts verified as "hits". Similarly, 
forty-eight hour wind speed forecasts verified 78.5; from isentropic 
chart data and 39% from 500 millibar chart data. So, percentagewise, 
forecasting wind vectors by CAV trajectories on isentropic charts 
was superior to forecasting wind vectors by CAV trajectories on 50C 
millibar charts during the period of this study. 

Figures 1 through 4 are histograms of error between observed 
and forecast wind speed and direction. Apeigaiie standard statistical 
methods, normal curves were fitted to the data. the test was 
applied to determine if the distributions of error were likely to be 


from a normal population. Only the twenty-four and forty-eight hour 


(13) 
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errors were compared and tested as the seventy-two hour isentropic 
sample was too small. In each case tested, the distribution was 
shown to be reasonably from a normal population. In cases of data 
from isentropic charts, there was little doubt that the distributions 
were from normal populations. But in cases from the 500 millibar 
charts, the values of Y~ were close to the critical values and, 
therefore, there is some doubt that they may be from a normal popu- 
lation. However, since the values of Y" did fall below the critical 
values, we will treat such data as being from normal populations. 

The Myr test [5, pp 152-154 | is a statistical test for comparing 
variances (s*) of two samples of data to determine if the two samples 
can reasonably be from the same normal population. The basic require- 
ment of the "F" test is that the data samples must each be from normal 
populations, although not necessarily the same normal population. This 
is why it was necessary to show normality of the data samples above. 

To apply the "F" test, we first postulate that there is no significant 
difference between samples tested, i.e., 5° = 3. . This is the null 
hypothesis. If the value of F computed falls below a critical value, 
the null hypothesis is not denied and there is no statistical difference 
between the two samples. If the value of I computed is higher than the 
critical value, the null hypothesis is denied, and there is a statisti- 
cally significant difference between the two samples and they are not 
reasonably from the same parent normal population. It is this last 
condition we wish to obtain here. The results of applying the "rr" 


test are as follows: 
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Twenty-four liour Wind Direction bLrror. The variance of 


data from the isentropic charts 5, = 1.24 =. The variance 
of data from the 500 millibar charts $.°= 6.37 . Then 
i= (Sag eye ere 
F. at a 95% level of belicf is 4.90. Therefore, the 
samples tested are significantly different and the null 
hypothesis is denied. 
Forty-eight iiour Wind Direction Error. The variance of 
date from the isentropic charts 5*= 3.15 - The 
variance of data fren the 500 millibar charts $= 30:72 . 
Then 
a SiS = 9.75. 
F. at a 95% level of belief is 3.75. This is highly 
Significant and the null hypothesis is certainly denied. 


Twenty-four Hour Wind Speed Error. The variance of data 





from the isentropic charts is 155.7. The variance of data 
from the 500 millibar charts is 413.1. ‘Then 

F= 2.65. 
Fe at a 953 level of belief is 5.19 and, therefore, the 
null hypothesis is not denied and there is no significant 


difference between the samples cf data. 


Forty-eight ilour Wind Speed lirror. The variance of data 


from the isentropic charts is 127.3. The variance of data 


from the 500 millibar charts is 515.5. Then 
F= 4.05, 


At a 905 level of belief the value of F,= 4.00, and at 
this level the difference between the samples becones 


Simmatieant. 


(15) 
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Significant results were obtained in the cases of wind direction 
error and in one case of wind speed error. In each of the cases 
where significance was obtained, an intuitive examination of the 
data shows tnat 1t is the data from trajectories constructed on 
isentropic charts that are better. 

It 1s gratifying to note tnat the greatest improvement was 
Shown for wind direction forecasting for a forty-eight hour period. 
In much operational forecasting, wind direction is usually deemed 
more important than wind speed for prognostic work. ‘rhe proper 
placing of troughs and ridges and associated weather patterns 
depends greatly on correct wind direction forecasts. Therefore, 
from the results obtained, the CAV trajectory method will show 
greatest improvement over the present system of taking CAV tra- 
jectories on constant pressure charts if used on isentropic charts 
to forecast wind direction, at least up to forty-eight hours. 

A recent study of CAV trajectories along a psuedo-600 millibar 
chart was published by Bruch [ 2]. Bruchts results on wind direction 
error are quite similar to results obtained herein for the 500 rullabar 
chart. Bruch did not verify wind speeds. 

It is realized that the isentropic chart leaves much to be desired 
&S an operational weather chart. At present the chart is of course 
laborious to plot since no isentropic data are transmitted on tele- 
type circuits; however, if such data were transmitted on teletype 
circuits, streamline analysis would take little longer to prepare 


than the present 500 millibar chart. Another difficulty is that tne 
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movement of the isentropic surface is not easy to predict. Since 
the CAY trajectory developec herein is a two dimensional projection 
of a three dimensional path, the height of the parcel in the future 
is in doubt. Of course, if the isentropic surface did not move, the 
future position of a parcel would be certain. But, this is not 
generally the case. A future study to correlate this study with the 
prognosis of isentropic surfaces could be to advantage. lNonetheless, 
the computation of CAV trajectories on an isentropic surface can be 
of value to prognosis of a nearby constant pressure surface. The 
positions of prognostic troughs and ridges of the isentropic flow 
ean be placed on a nearby constant pressure surface. 

From study of the CAV trajectories on the isentropic charts of 
this paper, it qualitatively appears that parcels imbedded in north- 
westerly flow descend as they move southeastward, whereas parcels 
inbedded in southwesterly flow rise as they move northeastward. 

In the final form of tne vorticity equation there is a diver- 
gence term that is assumed to be zero. ‘This assumption is not always 
fulfilled in the atriosphere. The effect of this divergence term may 
be important. For illustration, let the divergence term Vag" VY 
be constant in equation (9) and integrating with respect to¢ , we 


obtain 
~(%, Vv) ot 
(¥, +) - ce 


(17) 





a1 


Let {=O at the initial inflection point, where also i= 0° » then 
(f+) So f., 


~(V,,°V) dt 


> foe —f. (11) 


and 


As can be seen, positive values of the divergence term in a north- 
westerlv flow will cause a negative contribution to relative 
vorticity of a parcel causing it to go farther south than under a 
CAV trajectory. Positive values of the term in southwesterly flow 
will have a similar effect on the relative vorticity of a parcel 
causing it to not reach the maxinnm latitude indicated by a CAV 
trajectory. Negative values of the divergence term will cause 
positive contributions to the relative vorticity of a parcel. In 
northwesterly flow this will cause a parcel to curve more cycloni- 
cally than would be indicated by a CAV trajectory and to curve less 
cyclonically in southwesterly flow. Fultz (4, PP 32-109] mentions 
that in a lower level (10,000 foot level) the average type of 
deviation from CAV trajectories is such that could be caused by 
horizontal velocity divergence in northwesterly flow and convergence 
in southwesterly flow. Now, 

Veo’ = Ow 4+ ov — Ow pov 4}. Ow dz 4 Ov Dz 

OX oF \dy Of OZOY Oz OY 

Thus, the divergence term is composed of two parts, horizontal 


velocity divergence and a shear term. [Fleagle [3] and Panofsky [ 7] 


(18) 
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have given several values for horizontal velocity divergence, 
with a range of about 5 x 10-6 per second. The values of the 
shear term could be computed at initial inflection points for 
thirty-four of the CAV trajectories on isentropic charts of this 
study. In fourteen cases the term was less than 107? per s@¢ond. 
In eleven cases the value of the shear term was approxinately 10-°, 
and in nine cases, approximately -10~ per second. Usually, the 
Shear term was positive or negative if the flow was up or down the 
isentropic surface respectively. From this we see that in about 
60% of the cases studied, the shear term was of the same order of 
mascnitude as the horizontal velocity divergence. In these cases 
the shear term either canceled or doubled the horizontal velocity 
divergence term. Now, if the assimption of Woo: V =O were the 
most important assumption in the development of CAV trajectories, we 
might expect forecasting of future poSitions of parcels by CAV tra- 
jectories on isentropic surfaces to give, in some cases, greater 
error than from applying CAV trajectories to constant pressure sur- 
faces. This was not observed in the data collected in this study. 
This may suggest that the assumption of no horizontal divergence is 
not too restrictive, and that possibly one of the other assumptions 
made by Rossby is more restrictive, e.g., assumption of purely 
horizontal motion. 

If we take a large constant value for horizontal velocity 


divergence, say 6 x 1076 per second, take a like value for the 


(19) 
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shear term, asswne a tine period of twelve hours, and use these 
values in equation (11), we would find that the value of +. at, 
say, 45° is reduced by 33,5 or to the value of € at latitude 33°. 
But, this is an extreme case, and the effects of the divergence 
term are usually much less. However, since the assumption of 
Yeo 'V= © does lead to error, it seems that some future re- 
searcn could be profitably spent introducing quantitative values 
for the divergence term and computation of a trajectory incorpor- 
ating this correction to the CAV trajectory for use in more 


accurate forecasting. 
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Figure 2. 
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Figure 3. 
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